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The Bohigas-Giannoni-Schmit conjecture stating that the 
statistical spectral properties of systems which are chaotic in 
their classical limit coincide with random matrix theory is 
proved. For this purpose a new semiclassical field theory for 
individual chaotic systems is constructed in the framework of 
the non-linear cr-model. The low lying modes are shown to 
be associated with the Perron-Frobenius spectrum of the un- 
derlying irreversible classical dynamics. It is shown that the 
existence of a gap in the Perron-Frobenius spectrum results 
in a RMT behavior. Moreover, our formalism offers a way of 
calculating system specific corrections beyond RMT. 

The theory of random matrices 0] emerged from the 
need to characterize complex quantum systems in which 
knowledge of the Hamiltonian is minimal, e.g. complex 
nuclei. The basic hypothesis is that the Hamiltonian may 
be treated as one drawn from an ensemble of random ma- 
trices with appropriate symmetries. It has been proposed 
by invoking the complexity of systems which have many 
degrees of freedom with unknown interaction coupling 
among them. 

The study of the statistical quantum properties of sys- 
tems with small number of degrees of freedom, within the 
framework of random matrix theory (RMT), has devel- 
oped along two parallel lines. The first was by consid- 
ering an ensemble of random systems such as disordered 
metallic grains Randomness in this case is intro- 
duced on the level of the Hamiltonian itself, e.g. as a 
consequence of the unknown impurity configuration. In 
the second approach, RMT was used in order to under- 
stand the level statistics of non-stochastic systems which 
are chaotic in their classical limit such as the Sinai or the 
stadium billiards Here "randomness" is generated 
by the underlying deterministic classical dynamics itself. 
Nevertheless, it has been conjectured ||] that "spectrum 
fluctuations of quantal time-reversal invariant systems 
whose classical analogues are strongly chaotic have the 
Gaussian Orthogonal Ensembles pattern" . 

Despite being supported by extensive numerical stud- 
ies, the origin of the success of RMT as well as its domain 
of validity are still not completely resolved. In this letter 
we show that, in the semiclassical limit, this conjecture is 
indeed valid for systems with exponential decay of classi- 
cal correlation functions in time. Moreover, the formal- 
ism which we introduce below offers a way of calculating 
system specific corrections beyond RMT. 

So far, the main attempts to establish the relation- 
ship between non-stochastic chaotic systems and RMT, 
have been based on periodic orbit theory Gutzwiller's 



trace formula expresses the semiclassical density of states 
as an infinite sum over the classical periodic orbits of the 
system. However, the number of periodic orbits is ex- 
ponentially large and clearly contains information that 
is redundant from quantum mechanical point of view. 
This detailed information conceals the way of drawing a 
connection between the quantum behavior of chaotic sys- 
tems and RMT. Indeed, the success of the periodic orbit 
theory approach in reproducing RMT results |^ appears 
to be limited. 

Here we develop a new semiclassical approach in which 
the basic classical ingredients are not individual periodic 
orbits but global modes of the time evolution of the un- 
derlying classical system. It is possible to construct a 
field theory in which the effective action is associated 
with the classical flow in phase space. We argue that, 
the statistical quantum properties of the system are inti- 
mately related to the irreversible classical chaotic dynam- 
ics or, more precisely, to the Perron-Frobenius modes, in 
which a disturbance in the classical probability density 
of chaotic system relaxes into the ergodic distribution. 
These modes decay at different rates. This, in principle, 
enables a description of the system at levels of increasing 
complexity by incorporating higher and higher modes. 
The "zero mode" manifests the conservation of classi- 
cal probability and corresponds to a uniform distribution 
over the energy shell. We show that RMT coincides with 
a description which takes into account only this mode. 
The deviations from the universal RMT behavior emerge 
from the higher modes of the evolution of the classical 
system. 

Our approach is analogous to that of disordered sys- 
tems where the diffusion modes account for the classical 
relaxation. However, in the field theoretic description 
of disordered systems averaging is performed over an 
ensemble. Here, in contrast, in order to characterize in- 
dividual systems only energy averaging will be employed. 

The main result of this letter will be to establish the 
Bohigas-Giannoni-Schmit conjecture. To do so we will 
first show that quantum statistical correlators are de- 
scribed by a functional non-linear cr-model. Its low lying 
modes are identified with the Perron-Frobenius cigcn- 
modes of the underlying classical dynamics. Then we 
argue that, provided classical correlation functions decay 
exponentially in time, there is an energy domain where 
the zero mode contribution governs the behavior. This 
follows from the fact that in such systems the Perron- 
Frobenius spectrum has a gap. Finally we establish the 
relation to RMT by identifying it with the zero mode of 
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the constructed field theory. 

Let H define a quantum Hamiltonian whose classical 
counterpart, H{x), is chaotic. Here we use the notation 
X = (q, p) to denote a vector defined in a 2 x d dimen- 
sional phase space. We restrict attention to a system 
in which all classical trajectories are unstable, and there 
are no islands of regular motion in phase space. At suf- 
ficiently high energy Eq, the mean spacing A between 
adjacent energy levels can be assumed constant and ex- 
pressed through the Weyl formula. 
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Henceforth energy and time will be measured in units of 
A and the Heisenberg time h/ A respectively. 

The energy averaging is taken over an energy band 
containing a large number of levels N such that 1 <^ 
N <^ eo ~ Eo/ A. For simplicity we choose to work with 
Gaussian averages (e = E/A): 
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The basic quantity calculated within the field theory 
approach is the generating function {Z{J))^g. Any n- 
point correlation function can be obtained by taking its 
derivatives with respect to the various components of the 
source J. However, to keep the discussion simple, we will 
restrict attention to systems belonging to the unitary en- 
semble (i.e. those with broken T-invariance) , and focus 
on two-point correlation functions such as the correlator 
of density of states, R{s) — (p(e + s)/9(e))eo — 1, where 
p{e) = Tr(5(e — H). As long as s is sufficiently small 
compared to the bandwidth iV, the final results are inde- 
pendent of the particular form of energy averaging. 

Following the usual approach [Zj, the generating func- 
tion can be expressed as a field integral. Introduc- 
ing the four component supervector field \E'^(q) — 
("0^1 X^i V'"^X^)j where "0 and x denote commuting and 
anti-commuting components respectively, and the super- 
script A {R) denotes the advanced (retarded) compo- 
nents ifnl , we can define 



Z(J) = Jv-^e 
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where the action is given by 



5[*,e] =z / A'l'(q) 
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*(q). (4) 



Here s+ = s + iO, while A — diag(l, 1, — 1, — 1) and 
k = diag(l, — 1, 1, — 1) break the symmetry between re- 
tarded/advanced field components and supersymmetry 
respectively, and ^ = \1/^A. The use of fermionic as 
well as bosonic components ensures the normalization 
Z(0) = 1. To evaluate the two point density correlator 
one can choose J to be constant, then 
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Energy averaging (^) of Z{J) generates a quartic in- 
teraction of the 5* fields: 5(5', e) 5(5*, eo) -I- 5j, where 



^^i = ^ ( / rf'9*(q)*(q) 



(6) 



In contrast to impurity averaging 0, energy averaging 
induces a non-local interaction of the ^ fields. This in- 
teraction term can be decoupled by means of Hubbard- 
Stratonovich transformation, with the introduction of 
4x4 supermatrix fields (3(q, q') which depend on two 
coordinates, 
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Here STrq denotes the trace operation for supermatri- 
ces, while the subscript q implies a further extension of 
the trace to include integration over all spatial variables, 
e.g. STrqQ2 ^ / d'*gd'^(7'STrQ(q, q')g(q', q). Integrat- 
ing over ^! we obtain: 



VQ exp 
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g-\Q)^eo-Y^-H -JkA-NQ. (9) 



The derivation of Eq. (|^) involved no approximation. 
Further progress is possible only within a saddle-point 
approximation, which is accurate to order Vary- 
ing the total action with respect to Q one obtains the 
saddle-point equation, 



Qog'\Qo)=N 
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where Qq and GiQo) are operators. 

To understand the structure of the saddle-point man- 
ifold, in the scmiclassical limit, it is useful to employ 
the Wigner representation of operators. Given an op- 
erator O as a set of matrix elements 0(qi,q2) between 
two position states at qi and q2, its Wigner represen- 
tation is a function of the phase space variables x de- 
fined by ©(x) = / d''g'exp(ipq7;i)0(q-Hq72,q-q72). 
We will use the fact that, in the semiclassical limit, the 
Wigner transform of a product of operators is equal the 
product of the Wigner transformed operators, namely, 
(C'i02)(x) C'i(x)02(x), where e'i,2(x) are smooth 
slowly varying functions on quantum scale. 

Treating s and J as small compared to the bandwidth 
and introducing the phase space variables, X|| on 
the energy shell and x± — H perpendicular to the en- 
ergy shell, the solution of Eq. ( p^ can be written in the 
Wigner representation as 
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This solution is, however, not unique. In fact, the saddle- 
point solutions form a degenerate manifold in super- 
space associated with the group of pseudo-unitary ro- 
tations. Both the integration measure and the action in 
Eq. arc invariant under the group of transformations 
Q(x) U^^{x±)Q{x)U{x±), where U{x±) belongs to 
the pseudo-unitary super-group t/(l, 1/2) Thus, any 
matrix of the form Q{x±) = U^^{x±)Qo{x±)U{x±) is a 
solution of the saddle-point equation (||) . 

When integrating over the fluctuations 5Q, near the 
saddle-point manifold, one has to take into account the 
anisotropy of the dependence of the action on the fluc- 
tuations SQ. Those fluctuations on which the action de- 
pends strongly (massive modes) can be integrated out 
within a conventional saddle-point approximation. The 
remaining fluctuations describe the Goldstone modes of 
the system, and their integration must be performed ex- 
actly. The experience gained from studies of disordered 
systems suggests that these degrees of freedom can be 
parametrized by Q(x) — T^^{x)Qo(x±)T{x.). However, 
integration over the massive modes ||9| shows that, in the 
limit iV ^ 1, the only non- vanishing contribution comes 
from matrices T which are independent of the energy xj_ . 
The Goldstone modes can be therefore parametrized by 



g(x)=T-i(x||)Qo(2:±)r(x||), 



(12) 



to 



the coset space 



where T{x\\) belong 
C/(l,l/2)/C/(l/l)®C/(l/l) §. 

The derivation of the effective field theory can 
thus be obtained by: (i) Substituting Eq. (|l^) into 
Eq. (^), approximating commutators by Poisson brackets 
[C'i,C'2](x) i?i{Oi(x), 02(x)}, and replacing the trace 
by the phase space integral, Trq(C') — h~''' J (ixC'(x) 
(this is the entry point of the semiclassical analysis) (ii) 
Expanding the logarithm to first order in , J, and the 
Poisson bracket {iJ, T(x||)} (higher order terms of this 
expansion appear at order 1/7V and are neglected) (iii) 
Performing the x_l integration of the resulting action. 
The last step relies on the fact that within the energy 
band, where averaging takes place, the classical dynam- 
ics is independent of the energy. As a result we obtain 
the a-model: 
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Q = -^ J dx^T-H^\\)Qo{^±)T{^i\)^T-'AT, (15) 

and C is the dimensionlcss infinitesimal time evolution 
operator defined by the Poisson bracket 



£ ■ = n{H, ■}. 



(16) 



which relies on disorder averaging in the limit of vanish- 
ing disorder [^ . 

To make sense of the functional integral in Eq. ( p^ ) we 
must identify the low lying modes of the action. In the 
case of impurities, the low lying degrees of freedom cor- 
respond to the eigenvalues of the diffusion operator that 
form a discrete spectrum. In the present case it is tempt- 
ing to associate the low energy degrees of freedom with 
eigenmodes of the unitary (reversible) evolution operator 
e~^*. However this identification is incorrect. 

As with any functional integral there is a need to de- 
fine an appropriate regularization. For example the func- 
tional integral in Eq. (|l^) may be understood as the limit 
a ^ of a product of definite integrations over a dis- 
cretized space, where a denotes the discretization cell 
size. This admits to functions r(x||) which are smooth 
and square integrable. More generally, a regularization 
can be obtained using a truncation of an arbitrary com- 
plete basis. 

In seeking such a basis, the eigenfunctions of the clas- 
sical evolution operator C seem to be the natural choice. 
However, the intricate nature of chaotic classical evolu- 
tion cause these eigenfunctions to lie generally outside 
the Hilbert space. This can be seen from the follow- 
ing consideration: Chaotic dynamics of probability densi- 
ties involves contraction along stable manifolds, together 
with stretching along unstable ones. Thus, in the course 
of time, an initially non-uniform distribution turns into 
a singular function on the unstable manifold, which in 
turn covers the whole energy shell densely. Therefore the 
eigenfunctions of C, which require the infinite time limit, 
are not square integrable and their contribution to the 
functional integral can not be directly recovered by the 
discretization procedure involved in evaluating the func- 
tional integral. 

Thus, in choosing a convenient basis one has to take 
account of the regularization. Its primary effect is to 
truncate the contraction along the stable manifold and 
thereby render the classical evolution irreversible. It is 
the eigenfunctions of this regularized classical evolution 
operator that serve as a suitable basis for the quantum 
mechanical correlator. 

A natural way of introducing a regularization is to ex- 
press the functional integral in terms of a basis of eigen- 
functions of the evolution operator to which a small noise 
has been added such that a diffusion on the energy shell 
is allowed. The eigenfunctions of the resulting operator, 
which is no longer antihermitean, are rendered smooth 
along the stable manifold and square-integrable. When 
the strength of the noise is taken to zero — the spectrum 
of the resulting operator, known in the literature as the 
Perron-Frobenius operator, reflects intrinsic irreversible 
properties of the purely classical dynamics 
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approaches which recover this spectrum involved a use 



of symbolic dynamics |]12 
dynamics in phase space 



For a stochastic Hamiltonian, an action equivalent to 
Eq. (14) has been proposed recently using an argument 



course graining of the flow 
and analytic continuation 
methods ||lj]. We remark here that the physical spec- 
trum of the classical operator £ appears when it prop- 
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agates smooth probability densities. Thus the matrices 
T(x||) over which H, acts are understood to be smooth. 
In this respect, quantum mechanics can viewed as the 
natural framework for calculating irreversible properties 
of classically chaotic systems. 

Let {7n} be the set of the eigenvalues of the genera- 
tor, £, of the Perron-Frobenius operator. In ergodic sys- 
tems, the leading eigenvalue 70 = is non-degenerate, 
and manifests the conservation of probability density. 
Thus any initial density relaxes, in the course of time, to 
the state associated with 79. If, in addition, this relax- 
ation is exponential in time, then the Perron-Frobenius 
spectrum has a gap associated with the slowest decay 
rate. Thus, for the first nonzero eigenvalue 71 we have 
7^^ = Re(7i) > 0. This gap sets the ergodic time scale, 
Tc = 1/71, over which the classical dynamics relaxes to 
equilibrium. In the case of disordered metallic grains, it 
coincides with the Thouless time, while in ballistic sys- 
tems or billiards it is of order of the time of flight across 
the system. 

In the limit 5^7^ the dominant contribution to 
Eq. ( p^ ) comes from the ergodic classical distribution, 
the zero-mode LTq = 0. With this contribution alone, 
the functional integral (|l^) becomes definite, 

(Z( J))eo - y" dQo exp (-*|STr [(s+A + 2 JfcA) Qo] ) , 

(17) 

where Qo = TcT^ATq. This expression coincides with 
that for RMT and leads to the Wigner-Dyson statistics 

This result can be generalized to any n-point corre- 
lation function as well as to systems with T-invariance. 
One thus concludes that the quantum statistics of chaotic 
systems with exponential classical relaxation are de- 
scribed by RMT at energies smaller than 7^. 

RMT description is expected to hold even for certain 
chaotic systems where the Perron-Frobenius spectrum is 
gapless 1^, such as the stadium or the Sinai billiards 
where classical correlation functions decay algebraically 
in time [|6|. The resolvent 1/(2; — £), in this case, is 
expected to have cuts which reach the 3z axis. Never- 
theless, we expect the RMT description to hold when- 
ever the spectral weight of the resolvent inside the strip 
< 3?2; < 1 (which excludes the pole at the origin, how- 
ever) is much smaller than unity. 

The strength of the field theoretic approach is that it 
encompasses a range of energy scales which go well be- 
yond RMT. In particular, a similar procedure to that 
used in Ref. ||l^ can be employed to obtain -R(s) corre- 
sponding to the unitary ensemble 

1 3^ , , cos(27rs)„, , 

where T){s) = Yifj, ^fj-i^'^ + 7fi) "'^ the spectral deter- 
minant {Afj_ being regularization factors [|l8|) associated 



with the Perron-Frobenius spectrum. This confirms the 
conjecture made Ref. |l8|| . 

In conclusion, we have established the Bohigas- 
Giannoni-Schmit conjecture for chaotic systems in which 
classical correlation functions decay sufficiently fast. 
Moreover, the field theoretic approach allows the study 
of statistical characteristics of such systems on an en- 
ergy scale which is much wider than that in which RMT 
applies. These statistics are determined by the analytic 
properties of the classical resolvent operator l/{z — C). 
This theory, in principle, offers a systematic controlled 
way of investigating quantum corrections to the leading 
semiclassical description. 
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